Phonon-induced decoherence of the two-level quantum subsystem 
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■ Phonon-related decoherence effects in a quantum double-well two-level subsystem coupled to 

q 

C3 



O 

o 



a solid are studied theoretically by the example of deformation phonons. Expressions for the 
reduced density matrix at T = are derived beyond the Markovian approximation by means of 



explicit solution of the non-stationary Schrodinger equation for the interacting electron-phonon 



system at the initial stage of its evolution. It is shown that as long as the difference between 



the energies of the electron in the left and the right well greatly exceeds the energy of the 

m 

^ \ electron tunneling between the minima of the double-well potential, decoherence is primarily 

q ' 

(N| ■ due to dephasing processes. This case corresponds to a strongly asymmetric potential and 

O 

00 

spatially separated eigenfunctions localized in the vicinity of one or another potential minimum. 
In the opposite case of the symmetric potential, the decoherence stems from the relaxation 
processes, which may be either "resonant" (at relatively long times) or "nonresonant" (at short 
times), giving rise to qualitatively different temporal evolution of the electron state. The results 
obtained are discussed in the context of quantum information processing based on the quantum 
bits encoded in electron charge degrees of freedom. 

PACS numbers: 03.65.Yz, 63.20.kd, 73.21.La 
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I. INTRODUCTION 

Recent advances in fabrication of the quantum dot structures pQ, manipulation with the single 
atoms on a solid surface [2] , and atomically precise placement of single dopants in semiconduc- 
tors [Sill] make possible the construction of various solid-state architectures with predetermined 
characteristics. The state-of-the-art experimental techniques allow for a control of the quantum 
states of charge carriers in nanostructures [51 EE 13 El E] . A number of interesting phenomena 
have been observed and predicted, including Rabi oscillations [ID) , the entanglement between 
the states of interacting quantum dots [TT| . complete localization of the wave packet in one of 
the wells of a symmetric double- well potential perturbed by a monochromatic driving force [12J, 
the auxiliary- level-assisted electron transfer between the quantum dots [131 [H] , localization of 
two interacting electrons in a driven quantum dot molecule [15], coherent control of tunneling 
in a quantum dot molecule [16J, etc. 

The coupling of carriers to the surrounding crystal lattice results in entanglement between 
the carrier and lattice degrees of freedom and in the loss of coherence [17]. In experiment, 
decoherence of quantum states can lead to, e. g., decay of coherent optical polarization [10], 
damping of Rabi oscillations [18], errors in operations on the quantum bits (qubits) [ID] , etc. 
In theory, decoherence manifests itself in the temporal decay of the reduced density matrix 
elements. There exist different approaches to the description of decoherence effects in solids, 

see, e. g., Refs. [I7J [2D1 12D 1221 1231 123 • 

In this paper, we study the phonon-related decoherence of a two-level quantum subsystem 
within a solid by means of solution of the non-stationary Schrodinger equation for the interacting 
carrier-lattice system. We obtain an explicit expression for the state vector and find the reduced 
density matrix taking a trace over the phonon variables. Decoherence is shown to be primarily 
due to either dephasing or relaxation processes, depending on the specific structure of the 
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energy basis and eigenfunctions of the two-level subsystem. 

The paper is organized as follows. In Section II, we describe a model for the two-level subsys- 
tem coupled to a solid. In Section III, we present an approximate solution of the time-dependent 
Schrodinger equation for the interacting electron-reservoir system at zero temperature. An ex- 
pression for the reduced density matrix of the two-level subsystem is derived in Section IV, and 
several limiting cases are considered by the example of the double-dot structure. The results 
are discussed in Section V, and Section VI concludes the paper. 

II. MODEL FOR THE TWO-LEVEL SUBSYSTEM COUPLED TO A SOLID 

We consider an electron in a double-well potential formed in a solid. Such a subsystem 
corresponds to, e. g., the gate-engineered double-dot structure [23], two nearby donors beneath 
the semiconductor surface [26J, etc. We suppose that the two lowest states of the electron, |1) 
and 1 2), are well separated in energy from the excited states \k) with k > 3. Then the electron 
Hamiltonian can be written as 

#o = £i|l>(l|+£ 2 |2>(2| , (1) 

where E\ and E% are the eigenenergies of the stationary Schrodinger equation 

H \k) = E k \k) . (2) 

For the subsequent consideration, it is instructive to write the Hamiltonian ([T]) in the basis 
{\L), \R)} formed by the ground states of the electron in the left and the right well in the case 
that the wells are isolated from each other. We assume the wave functions (r|L) and (r\R) to 
be strongly localized in the vicinity of the corresponding potential minima. Then, neglecting 
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the overlap (L\R), we have 

H = E L \L)(L\ + E R \R)(R\ -j(\L)(R\ + \R)(L\) , (3) 

where El and E R are the energies of the states \L) and \R) localized in the vicinity of, respec- 
tively, the left and the right potential minimum, and A/2 is the energy of the electron tunneling 
between the two minima. 

The states |1) and |2) are related to the states \L) and \R) by the following expressions 

\1) = C_\L)+C+\R), \2) = C+\L)-C-\R) , (4) 

where 

1 



C± = — 
*2 



1± , El ~ Er (5) 
y/(E L - E R f + 

The relation between Ei }2 and El,r is 

E l + E rT ^(E l - E R f + A 2 
E i,2 = g ' 

In what follows, we shall consider the limiting cases of (i) the strongly asymmetric potential, i. 
e., E R — El >> A, so that we assume A = 0, an (ii) the symmetric potential, i. e., El = E R 
and A ^ 0. From Eqs. (jlj)-(jSD one has 

|1) = |L), \2) = \R) , E h2 = E L , R (7) 



in the case (i) and 

|l) = ' L > + ' fl >, |2) = ' L > - |fi > , Ei2= ^ + ^TA =g A (8) 

l / ^2 1 7 V2 ' 2 2 v ; 

in the case (ii). 

The phonon term in the Hamiltonian is (hereafter the Planck constant h — 1 if not stated 
explicitly) 

H Ph =i:^(bf b P+l) > ( 9 ) 



where ujp is the frequency of the phonon mode f3 = (q, A) with the wave vector q and polarization 
A, and bp (bp) is the creation (annihilation) operator of the phonon of the mode (3. The 
eigenstates ({n^}) of the stationary Schrodinger equation 



HMnp}) = E({n fi })\{n p }) 



(10) 



are defined by the set {rip} of the phonon numbers rip for every mode (3, the eigenenergies being 
equal to 

E({np}) = + 5) ■ (11) 



The electron-phonon interaction term is 



P 1 



(12) 



where p(q) = / dre iqr p(r) is the Fourier transform of the electron density operator p(r) = 
Hmn ^m( r )^n( r )\ m )( n \y an d * s the microscopic electron-phonon interaction matrix element, 
which can be expressed in terms of the deformation potential H and the density of the crystal 
p as (here we restrict ourselves to deformation phonons) 



A/3 = gS 



1/2 



2pujpV 



(13) 



with V being the normalizing volume. The Hamiltonian f[T2"j) can be written in the spin-boson 
form as 



Hi n t — &z 

P 



9pbp + 9*pbp 



(14) 



where a z = \L)(L\ — |-R)(-R| is the Pauli spin operator in the basis {\L), \R)}, 
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A, 



A L (q) - A R (q) 



(15) 



and A L , R (q) = f dre^\(r\L, R)f 



To be specific, in what follows we consider the double-dot system with the Gaussian-shaped 
electron wave functions (t\L, R) ~ exp(— |r — jj| 2 /2/ 2 ), where are the coordinates of the 
dot centers, and I is the effective dot size. Then from Eq. (1151) one has [X91, 



9(3 = ^exp I — I sin I — I , (16) 

where d = — tl (i.e., d is the interdot distance), and we choose the origin of the coordinates 
in between the dots. Note that the condition of vanishingly small overlap (L\R) implies that 
I « d. From Eqs. (TT5I) and (TTB1) we find the spectral density J{uj) that fully describes the 
effect of the phonon bath on the two-level electron subsystem 



^ 3 



sin 



UJ \LJ d J 



,2 



exp H^J' (i7) 



where we assumed the linear dispersion law cj q = sq with s being the sound velocity and used 
the notations u>d = s/d and lo\ = s/l from Ref. [28] (note that lo\ » u>d). 



III. STATE VECTOR OF THE INTERACTING SYSTEM 

The state vector \^(t)) of the interacting electron-phonon system satisfies the non-stationary 
Schrodinger equation 

<3f»=*l*W>, (is) 

where 

H = H + H P h + H int (19) 

is the full Hamiltonian. The state vector can be represented as a linear combination of products 
of the electron and phonon states, 

!*(*)> = E E C^ } (t)e^^^\k)\{n p }) . (20) 

k=l,2{n } 



Note that in the absence of the electron-phonon interaction, the coefficients C k ^ n& } do not 
depend on time. We consider the case of zero temperature, so that initially there are no 
phonons in the solid (np = for any mode (3), while the electron is in the superpositional state 
oll\L) + a R \R), where \o.l\ 2 + \&r\ 2 = 1, and hence the state vector at t — is 

(21) 



tf(0)>= (a L \L)+a R \R))\O ph ) 



where \0 p h) is the state without phonons, or, equivalently, 



Cx,o pfc (0) = oi L , C R)0ph (0) = a R , C fc ,{ n(3 ^o}(0) = . 



(22) 



is on 



[Relation between the coefficients Ci t { n „}(t) , C 2 ,{ ri/ ,}(t) and CL,{ n0 }(t), C R ^ nfj }{t) depends 
the relation between the states |1), |2) and \L), \R) which is different in the cases (i) and (ii), 
see Eqs. © and ©]. 

Substituting Eq. (|2"U|) into Eq. (TT5T) . we have a system of coupled differential equations for 
coefficients C k){ng} (t), 



~dt 



, ■ - e -iKfet-iS({n^})f 



E E Q, W} (t)(A;|({n /3 }|^|/)|{m /3 })e-^ i -^» f 

1=1,2 {m } 



(23) 



From Eq. ffT4l we find 



({n^KfclfT^IOKm/,}) 



W<^I0E 9p'\fmp> + 18 {mp },{ np }-ip, +^'v /7F ¥ 75 W},W}+V 
/3' L 



(24) 



where the designation {ng} Hp> means the set of the phonon numbers with one phonon of the 
mode (3' more/less than in the set {np}. With this expression, Eq. fl23|) becomes 



dt 



=1,2 



E WM>^ (M)t E c^^u^v^^ 



+G, { n, }+V (t)^ v /^' + le ^' < 



(25) 



where we took into account that E({np} ± 1^) = E({rip\) ± ujpi. From Eq. (|25|) we have 
t dC k , 0ph (t) 



E {fc|^|Z)e-^-^')^C, 1 ,(t)^e-^ (26) 



ut 1=1,2 P 



and 



z- 



Z=l,2 L /?' 



(27) 



We wish to consider either the initial stage of the system evolution or the case of "weak 
decoherence" , so that the state vector |\&(i)) does not differ much from |\l/(0)). This implies 
that all coefficients Ck,{ n/3 }(t) in \^(t)) except Ck >0ph (t) are small and decrease with the number 
of phonons J2p n p i n the state {np}, and hence, in the first approximation 

Mi)) « E C fc , 0p;i (t)e-^^V^)|0 M ) + E E^i,(*)^ £fc ^^^)|l^) , (28) 

fe=l,2 fc=l,2 p 

where Eq h = J2p^p/2 is the zero-point phonon energy. Neglecting the second (two-phonon) 
term in the right-hand side of Eq. ( 1271) . taking Ci^ ph (t) « Cy) fe (0), and integrating Eq. ( 1271) . 
we have 

e -i{Ei-E k -u g )t _ i 

<W*) « # E a^(0)<A;|^|/>-p p • (29) 

1=1,2 n,i- h, k -up 

Substituting this expression into Eq. ( f26l) and evaluating the integral, we have 

c fc) o, fc (*)«Cfc,o ph (o)+ E E c^(o)<*|^|/><^,N 

!=l,2ra=l,2 

1 / e -i{ E m-E k )t _ i e -i(Ei-E k +u )t _ jx 
^^ |2 F P^Tl - F F F TP + J" (30) 

Now let us analize the expressions for Ck,a ph {t) and C^i^i) in the limiting cases (i) and (ii) 
mentioned in Sec. II. 



(i) Strongly asymmetric double-well potential (El ^ Er,A = 0) 
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In this case, we have from Eqs. (J7J) 



[1\& Z \1) = 1, (2|a 2 |2> = -1 , (l\a z \2) = (2\a z \l) = 



(31) 



so that Eqs. ([29]) and ([30]) become 



C k>lg (t) « (-l) fe C Mph (0)^ 



(32) 



c k)0ph (t) * C k)Oph (0) 



1 + E + ^ - 1 



(33) 



It follows from Eq. (1331) that the quantitative conditions for our approximation Ch,o ph (t) 
C k,o ph (°) are: 



tfup Jo 



A J M 



and 



t << t r 



\m 



7 j f- j 







ClLd— — 



(34) 



(35) 



For the spectral function given by Eq. f[T7]) one has (the second term in square brackets of Eq. 
ffTT]) can be neglected since uj\ » Ud): 



^2 



A 



2 . ' 



8ir 2 ps 3 Ph h 2 ut 



^ps 2 l 3 h 1 
t » 87rV27r^ 



(36) 
(37) 



For the parameters of GaAs (S 7 eV, s = 5.1 • 10 5 cm/s, p = 5.3 g/cm 3 ) we have A ~ 3 • 10 -4 
and to ~ 10~ 8 s at / = 25 nm. Note that A decreases with I, while to increases, so that our 
approach works better (and in a broader time interval) for relatively large quantum dots. For 
example, A w 2 ■ 1CT 5 and to ~ 7 ■ 1CT 7 s at I — 100 nm (approximate dot size in Ref.[6]). 



(ii) Symmetric double-well potential (El = Er,A ^ 0) 
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Now it follows from Eqs. © that 



(l|<x 2 |l> = {2\a z \2) = , (1|^|2) = (2\o z \l) = 1 



(38) 



and from Eqs. ( 1291) and (1301) we have 



C 2 , l0 (t) * -C lj0ph (0)gp 



e i(u g -A)t _ 
Up — A 
e i(u>g+A)t _ J 



Ci, 0p ,(t) « Ci,o*(0) 
Ca^Ct) ~ C 2j0ph (0) 



1 + ?(^ + A)^ 
In I 2 

i + y^M 

. + ^{^-ay 



UJ/3 + A 

+ A)t + e - l ^ +A)t - 1 
- A)t + e - t{ ^- A)t - l) 



(39) 

(40) 
(41) 

(42) 



At A << uji (i.e., for well separated quantum dots), the conditions that coefficients Ckfl ph (t) 
do not differ much from their initial values Ck t o ph {0) are the same as in the case (i), see Eqs. 
(El and (ESI). 



IV. DENSITY MATRIX 

Having found the state vector \^(t)) of the interacting electron-phonon system, we can 
calculate the reduced density matrix of the two-level electron subsystem by tracing out the 
phonon variables. From Eq. ([20]) we have 



pit) = Tr {ng} m))m)\ = E E E e- i(Efc - Ei) *C &){n3} (t)C* { ^ } (t)|A;)(i| . (43) 

{np}k=l,2 1=1,2 

According to approximations made in Sec. Ill, see Eq. ( l28l) . we retain in the sum over {np} 
the zero- and one-phonon terms only, so that the matrix elements of p(t) in the energy basis 
read 



Pki(t) ~ e 



-i{E k -Et)t 



c k , 0ph (t)ci 0ph (t) + £ c kjl0 (t)c^ (t) 
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(44) 



Following the consideration in Sec. Ill, below we calculate Pki{t) in two limiting cases (i) and 
(ii). 



(i) Strongly asymmetric double-well potential (El Er,A = 0) 

Substituting expressions fl32l) and fl33|) into Eq. f|44|) . and neglecting the terms containing 
\gp\ A (i.e., the terms of the order of A 2 , (t/t ) 2 , and A(t/t )), we have 

p[t) ~ { p 21 (0) [1 - B 2 (t)] e -«M0' p 22 (0) ) ' l45j 

where p fci (0) = C fcj0jA (0)C^(0) and 

(Here we use the notations of Refs. pJ2, EH, EH]). The value of B 2 {t) equals to zero at t — 
and increases with t as B 2 {t) « 4A(^t) 2 at t << c^f 1 (the value of cjf 1 is about 2 • 1CT 11 s at 
I = 100 nm) up to B 2 (t) m 4A at t » uf 1 (this is consistent with the range of validity of our 
approximation, t « t , since t ~ 1/Aui » ouf 1 ). 

Since in this case |1) = \L) and |2) = \R), we arrive at the conclusion that electron-phonon 
interaction has no effect on the diagonal matrix elements pn and Prr, while the absolute values 
of the non-diagonal matrix elements pm and Prl decrease with time. So, there is no relaxation 
in the system (i.e., occupations of both dots do not change during the system evolution), and 
decoherence emerges as pure dephasing, in accordance with the results of Ref. [19] and "semi- 
classical" approach to electron-phonon interaction in the two- level systems [3D] ■ The reason for 
this is that in the case of strongly asymmetric double-well potential, the electron Hamiltonian 
ifo commutes with the interaction Hamiltonian Hi nt , see Eqs. ([1]) and (JHj). Contrary to the 
results obtained in Refs. [2TJ [29], our expression (145]) for p(t) contains the factor 1 — B 2 (t) 
instead of e~ B2 ^ . This is a consequence of our approximations. Indeed, B 2 {t) < 4A << 1, so 
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that e B ® wl - B 2 (t). For times t >> ^ 1 we have B 2 {t) rs 4A, i.e., the dephasing process 
does not depend on time. 

(ii) Symmetric double-well potential (El = Er,A ^ 0) 

Substituting expressions fl39|) . f|40|) . (jUI) . f|42l) into Eq. (jUJ, and neglecting, as above, the 
terms of the order of \g/3\ A , we obtain the density matrix elements in the energy basis: 

Pn(t) = 1 - pn(t) « Pn(0)fl -4F+(t)) + 4p 22 (0)F_(t) , (47) 



Pn(t) = p* 21 (t) a Pl2 (0)e iAt 
l^l 2 

; (w/j-A) 



, (^ + A) 2 



1 + ^ '^ |2 ( + A)t + e -^+ A )< - l 



+2p 2 i(0)E 



2 



COs(At) — COs(u!/3t) 



(48 



where 



'^-Sd^-"^)- (49) 

First we analyze the diagonal matrix elements Pkk{t)- For A << uj\ (for example, in ex- 
periment [6] the typical value of A is 10 10 s _1 , while Wj « 5 • 10 10 s" 1 , and A quickly 
decreases with the interdot distance d) we have F + {t) « B 2 (t)/8, see Eq. ( f46l) . As for 
F_(t) = J °° dujJ{uj) sin 2 [(c<j — A)t/2]/(w — A) 2 , one can distinguish two different contribu- 
tions to this term, one being from the "resonant component", i. e., from the 5-function-like 
peak of sin 2 [(a; — A)t/2]/(u — A) 2 as a function of u at uj — A, with a height of t 2 /4 and a 
width of ~ 1/t, and the other from the remaining "nonresonant background" of the phonon 
spectrum. Making use of the expression 

sin 2 (et) 



irte 2 
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S(e) , (50) 



the former can be estimated as F_(t) m Tt/A, where 



POO 

F = 2vr \9ib\ 2 S(uJi3 - A) = 2tt J duoJ(uj)5(uj - A) = 2vrJ(A) 



(51) 



is the rate of the phonon- induced electron transitions |2) — > |1) (relaxation rate) calculated by 
the Fermi golden rule [31], and the latter as F^\t) « B 2 (t)/8, so that « n/4 + 5 2 (t)/8. 

Note the different physics behind the two contributions to F_(t). While the term F!i\t) reflects 

(2) 

the energy conservation, hup = A, for a transition |2) — - > |1), the term F_ (t) arises because of 
the violation of the energy conservation at short times [31] . 
Since pn(0) = 1 — P22(0), one has from Eq. (11711 : 



# 2 (t) 



i - rt - 5 2 (t) 



-p 22 (0)e- r *- B2 W 



(52) 



where we took into account that B 2 (t) << 1 at any t and Ft << 1 at sufficiently short times (in 
fact, the standard first-order perturbation theory for calculation of the transition probability 
is valid at Ft << 1, see [31]). This expression for pn(t) differs from the commonly used one 
[20] by the presence of B 2 (t) terms which are responsible for the " non- resonant relaxation" and 
are time-independent at t » tof . The "resonant" term Ft prevails at very short times (note, 
however, that approximation (1501) which leads to the linear dependence of F_ on t may not 
hold in this time domain) and, more importantly, at t > t « A/J(A). For ujd « A << uj\ 
from Eqs. (H7J) and fl3SD one has t w (uJAfcuf 1 w 2.5 • 10~ 9 s at u; ; /A = 5 and I = 100 nm. 
So, the exponential changes in p n and p 2 2 go into play at t ~ t, while at tu^ 1 « t « t the 
diagonal matrix elements are, to the first approximation, time-independent (but different from 
their initial values), being determined by the value of A. 

As for the non-diagonal matrix elements, from Eq. (|4*8l) at A << u)\ we have 



Pu(t) ~ Pi2(0)e 



iAt 



l-2F + (t)-2F_(t) 



+ 2p 2 i(0) 



Acos(At) - A + 



S 2 (i) 
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Pn(0)e iAt 1 



B 2 



(t) Tt 



+ 2p 21 (0) Acos(Ai) - A + 



B 2 (t)] 



(53) 



2 



2 
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Just like pn(t), along with the usual "resonant" term Tt, this expression contains the "non- 
resonant" terms B 2 (t) and A. At cof 1 « t « t these terms are greater than the "resonant" 
one and thus govern the evolution of non-diagonal matrix elements. It is interesting, however, 
that upon going from the matrix elements pu in the energy basis {|1), |2)} to the matrix 
elements in the basis {\L),\R)} of localized electron states, the "non- resonant" contributions 
at t >> cof 1 can cancel out. For example, if electron initially occupies the left well of the 
double-well potential, |^(0)) = \L), i.e., Pll(O) = 1, Plr(0) = Pkl(O) = Prr(0) = and hence 
Pn(0) = Pi2(0) = P2i(0) = P22(0) = 1/2, then the probability to find an electron in the right 
well at time t is 



This simple expression describes the damped electron oscillations between the left and right 
wells and agrees with that obtained by Wu et al. [28J. 




Solid-state systems are of great interest in searching for a scalable quantum computer tech- 
nology, see Refs. [32J [33] and references therein. In particular, spatially separated orbital 
states of an electron in a pair of tunnel-coupled quantum dots [13], [3U [35] (or in a singly ionized 
pair of phosphorous donors in silicon [26], etc.) can be used as logical states of a quantum 
bit (qubit), the logical |0) (|1)) being associated with the state \L) (\R)) localized in the left 
(right) double- well potential minimum. These so-called charge qubits can be manipulated, e.g., 




(54) 




V. DISCUSSION 
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by applying adiabatically switched gate voltages [6j [7J [H [9j [26] or laser pulses (both resonant 
and off-resonant) [T31 ESI EZl ESI EH] to the system. 

The phase gate is realized if the energies El and E R of the states \L) and \R) differ from 
each other, while electron tunneling between these states is suppressed. Then the qubit vector 
evolves as = C L (t)\L) + C R (t)\R) = C L (0)e- iELt \L) + C R (0) e - iERt \R) , the absolute values 

of Cx(t) and C R (t) remaining unchanged, and the relative phase i(-Ei? — E£)t varying linearly 
in t. This corresponds to the case (i) of strongly asymmetric double- well potential considered 
above. As follows from the expression (|45p for the density matrix, the deformation-phonon- 
induced qubit decoherence during the phase operation is entirely due to dephasing processes 
and is quantified by the value of B 2 (t), see Eq. fj46l) . At t >> ujf 1 , the value of B 2 {t) becomes 
time-independent and equals to a constant 4A, which depends on the material parameters (p, 
s, H) and the quantum dot size I. For example, this constant is smaller than 10~ 2 in the GaAs 
based quantum dots with I > 10 nm, see Eq. (I3"6"j) . 

For the amplitude gates, the weights of \L) and \R) states in the qubit state |\&) change 
with time (the combination of the phase and amplitude gates allows for an arbitrary rotation 
of the qubit vector on the Bloch sphere). For example, \Cr{T)\ = |Cx(0)| and |Cl(T)| = 
|Cr(0)| at operation time T for the quantum NOT. In the absence of decoherence, this gate 
is implemented at El = E R (i.e., in the case (ii) of symmetric double-well potential) in time 
T = irk/ A m 2 • 10~ 10 s at A px 10 /xeV, see Eq. (154^) . Decoherence results in the damping 
of coherent electron oscillations between the dots. At the very early stage of qubit evolution, 
t < uf 1 = l/s (e.g., at t < 2 ■ 10~ n s for the GaAs dot size / = 100 nm), both "resonant" and 
"nonresonant" relaxation processes contribute to the decoherence, see Eq. ( 153)) . Contrary to 
the phase gate, at t > uf 1 the decoherence is primarily due to usual "resonant" relaxation [40], 
and the gate fidelity decreases with time as 1 — 0.5 exp(— Tt/2), see Eq. f[54l . The value of T is 
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extremely sensitive to the system parameters such as the dot size, the interdot distance, etc., 
see Eqs. (JT7J) and (15 ip . For example, at A >> Ud, the decoherence rate T oc A 3 exp(— A 2 /2uf) 
first increases with A up to the maximum value T max and next decreases rapidly. The value of 
r mQ x is about 5 • 10 9 s _1 in the GaAs dots with I = 10 nm and decreases with I as T max oc /~ 3 . 
At A << the decoherence rate is very small, V oc A 5 . 



VI. CONCLUSIONS 

In this paper, we analized the effect of deformation phonons at the initial stage of coher- 
ent electron dynamics in the double-dot structure by the examples of symmetric and strongly 
asymmetric double-well potential. We have explicitly shown that the phonon- induced decoher- 
ence can be due to both dephasing and relaxation ("resonant" and " nonresonant" ) processes, 
the decoherence rate being determined by the material and double-dot parameters. Making 
use of the appropriate spectral function J(u), the results obtained can be applied to describe 
the decoherence due to electron coupling with piezoelectric phonons in the double-dot sys- 
tem [19], [281 HD H21 HIH Hi] and with acoustic phonons in the double-donor Si-based structure 
[T9| 1301 05]. Generalization to the case of nonzero temperature (HI H21 [J3l 05] is straightfor- 
ward. To study the non-Markovian electron dynamics in more detail, it would be interesting 
to extend the consideration to the longer evolution times through account for A-phonon states 
with A > 2. 
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